Charge noise at Cooper-pair resonances 
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We analyze the charge dynamics of a superconducting single-electron transistor (SSET) in the 
regime where charge transport occurs via Cooper-pair resonances. Using an approximate description 
of the system Hamiltonian, in terms of a series of resonant doublets, we derive a Born-Markov master 
equation describing the dynamics of the SSET. The average current displays sharp peaks at the 
Cooper-pair resonances and we find that the charge noise spectrum has a characteristic structure 
which consists of a series of asymmetric triplets of peaks. The strongest feature in the charge noise 
spectrum is the triplet of peaks centered at zero frequency which has a peak spacing equal to the 
level separation within the doublets and is similar to the triplet in the spectrum of a driven, damped, 
two-level system. We also explore the back-action that the SSET charge noise would have on an 
oscillator coupled to the island charge, measurement of which provides a way of probing the charge 
noise spectrum. 

PACS numbers: 74.50.-|-r, 74.78.Na, 42.50.Lc 
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I. INTRODUCTION 

Mesoscopic superconducting circuits in which there is 
an interplay between Josephson tunnehng and electro- 
static charging effects display a wide range of interest- 
ing behaviors. The quantum coherence and non-linearity 
generated by Josephson junctions make superconducting 
circuits obvious candidates for qubitsi*^ However, the 
usefulness of qubits depends strongly on the extent to 
which the level of dissipation and the attendant deco- 
herence can be minimized.^ In contrast, there are many 
other applications of superconducting circuits, as mea- 
suring devices^^ or as coolers,^ where irreversible trans- 
port is necessary and hence dissipation plays an essential 
positive role. 

The superconducting single electron transistor^ 
(SSET) is an example of a superconducting device with 
applications where both coherence and dissipation are 
important. The SSET consists of a small superconduct- 
ing island linked to superconducting leads by Josephson 
junctions. A voltage gate coupled to the island can be 
used to control the flow of charge. The strong depen- 
dence of the charge transport through the SSET island 
on the properties of the gate makes it an ideal sensing 
device; it can act as either an electrometer^^ or, (when 
the gate is mechanically compliant) as a displacement 
detector. '^^^ Sensing typically involves a coupling between 
the SSET island charge and the degree of freedom be- 
ing measured, hence the charge noise spectrumiii^iii 
determines the back-action of the SSET on the mea- 
sured system (within the linear response regime). Al- 
though the back-action is a nuisance in the context of 
measurement, it can be used to manipulate the state of 
the measured system )i^ii^ii^ one recent experiment used 
the back-action from a SSET to cool a nanomechanical 



resonator,^ whilst another demonstrated the production 
of laser-like states of self-sustained oscillation in an elec- 
trical resonator. ^^ From another perspective the back- 
action of the SSET can be seen as providing an efficient 
means of probing its quantum noise properties ^li^ 

The SSET supports a wide range of different current 
carrying processes depending on the choice of operating 
point^~— When a large voltage is applied (V > 4A/e, 
where A is the superconducting gap), it becomes ener- 
getically favorable to break up the Cooper-pairs and the 
transport is dominated by inelastic (incoherent) tunnel- 
ing of quasiparticles at both junctions. At lower voltages 
the transport involves Josephson tunneling of Cooper- 
pairs, though dissipation is still required to generate a 
dc current. In the regime where V ~ A/e, current 
resonances known as the Josephson quasiparticl o^^i^^i^^ 
(JQP) and double Josephson quasiparticle^^^^ (DJQP) 
cycles combine coherent Cooper-pair tunneling and dis- 
sipative quasiparticle tunneling. The JQP and DJQP 
cycles have attracted much recent attention as the res- 
onant and (partly) coherent nature of the transport 
leads to measurement sensitivities approaching the quan- 
tum limit, as well as a range of interesting back-action 
effects.4i^aaa2-i5^^ 

In this article we focus instead on Cooper-pair 
resonancesii — i^i2^ — which occur at even lower volt- 
ages than the JQP and DJQP cycles and involve the 
coherent transfer of one or more Cooper-pairs across the 
two SSET junctions. At these lower voltages quasipar- 
ticles are almost completely absent and dissipation and 
decoherence are dominated by the electromagnetic envi- 
ronment of the SSET. Cooper-pair resonances are known 
to give rise to sharp features in the current voltage char- 
acteristics of the SSET,ii--iai^i26-i28 though much less is 
known about their quantum noise properties. 



We develop a simple model for the SSET system, valid 
at low voltages and for embedding circuits with resis- 
tances much less than the resistance quantum, Rq = 
h/Ae^. This model, which has important similarities 
with the standard description of resonance fluorescence 
of quantum optics, divides the energy levels of the 
system into pairs of resonant levels with the spacing 
within a doublet much less than the spacing between 
doublets^i22 We derive the master equation of the SSET 
within the Born-Markov approximations and then calcu- 
late both the average current and the charge noise spec- 
trum in the vicinity of the Cooper-pair resonances. The 
dominant feature in the charge noise spectrum is a triplet 
of peaks centered at zero frequency that is a character- 
istic of a driven two-level system, ^^'^^ additional weaker 
triplets of peaks also occur at much higher frequencies. 
We also investigate the back-action that the Cooper-pair 
resonances would have on an oscillator coupled to the 
SSET island. The asymmetry in the triplet of peaks cen- 
tered on zero frequency means that the resonances could 
be used to cool a mechanical resonator with frequency in 
the MHz range though not all the way down to its ground 
state. 

The organization of this paper is as follows. In Sec. 
im we describe our model of the SSET at low voltages, 
discussing the origin of the Cooper-pair resonances and 
the way in which the electromagnetic environment cou- 
ples to the charge passing through the transistor. Then 
in Sec. IIIII we show how the Hamiltonian of the SSET 
close to a Cooper-pair resonance can be transformed to 
give a description in terms of a set of doublet states. 
We then use standard approximations to derive a master 
equation that describes the SSET charge variables in the 
presence of dissipation due to the electromagnetic envi- 
ronment. Next, in Sec. lIVl we use the master equation to 
derive the average current at the resonances. The charge 
noise spectrum of the SSET is calculated in Sec. |Vl and 
the back-action is discussed. Finally we give our con- 
clusions in Sec. IVII The appendix contains additional 
details about certain aspects of the calculations. 



II. MODEL SYSTEM 

The SSET is shown schematically in Fig. [H For sim- 
plicity, we assume that the drain-source bias, V ^ is ap- 
plied symmetrically and take the junctions to have equal 
capacitances, Cj, and Josephson energies, Ej. A volt- 
age, Vg^ is applied to the gate which has capacitance, C^, 
(assumed to be much less than Cj). The electromagnetic 
environment of the SSET is modeled by the impedance 
Zq. We assume that the charging energy of the island, 
Ec = 4e^/2Cs, with C^ = 2Cj ^ Cg, is the dominant 
energy scale in the system, so that Ec ^ k^T^ Ej^ where 
T is the temperature. In this limit the SSET is best de- 
scribed in terms of charge states, using the two quantum 
numbers, n = Nl — Nr, the excess number of Cooper- 
pairs on the island, and k = {Nl + Nb)/2^ the average 
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FIG. 1: Circuit diagram of the SSET. The SSET island is 
linked to the leads by the Josephson junctions, Jl(r), and is 
coupled to the voltage gate by the capacitance Cg. The bias 
voltage, V, and impedances in the circuit, Zo, are taken to be 
distributed symmetrically. 



number of Cooper-pairs which have traveled through the 
system, with A/'j^(i?) the number of Cooper-pairs which 
have passed through the left (right )-hand junction of the 
SSET. — i^i2S The large charging energy means that only 
the two states with values of n closest to the gate induced 
polarization charge, Ug = CgVg/2e^ have an appreciable 
chance of being occupied. 

The full Hamiltonian of the system can be written as. 



H = Hs + Hi, 



H, 



(1) 



where i^env is the Hamiltonian of the electromagnetic 
environment and the SSET Hamiltonian, Hs^ consists of 
two parts Hs = Hch + Hj. The charging Hamiltonian of 
the island is 

^<=h= ^ ^[^c(n-n<,)2-2eT/A:]|n,fc){n,fc|, (2) 

n=0,l k 

where we have taken < n^ < 1, so only the states |0, /c) 
and |1, /c) need to be considered. 

The states {|n, /c)} separate into two ladders: 

{|0,fc),|l,fc + 1/2)} and {|0,/c + V2),|l,^)}, where k is 
now an integer. Josephson couphng links together adja- 
cent members of the same ladder, but does not connect 
states from different ladders. Assuming that quasiparti- 
cle tunneling (which can link the ladders) is negligible, 
we need only consider one of the sets of states. Choosing 
the ladder {|0, /c) , |1, A: + 1/2)}, the Josephson coupling 
between states is given by^i^S 

Hj = - J ^(|0, k) (1, k + V2I + |0, k) (1, k - V2I + h.c), 

k 

(3) 
where J = Ej/2. 

We assume that the dominant source of dissipation and 
decoherence of the SSET is the electromagnetic environ- 
ment, modeled by the impedances in the leads connecting 
it to the voltage sources. The impedances lead to fluctu- 
ations in the drain-source voltage, JV, which couple to 
the system operator, /c, and give rise to the interaction 
Hamiltonian, 
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FIG. 2: (Color online) Energy levels of ifch, (a) with V = ^1^ 
at the p = 1 resonance and (b) with V > VreJ for p = 0. The 
dashed lines enclose the doublets: the energy levels which are 
almost degenerate near the given resonance. 



The effects of 5V on the SSET are determined by the 
spectrum of the voltage fluctuations which takes the 
form^ 



Sv{oj) 



4e^ / {6V{t)5V{0))e'"^ dt, 

^ — CO 



]_ Q—huj/k-BT 



Re[Zj 



(5) 



where Zt = {Z^^ — iujC^j j Cy)~^ is the total effective 
impedance seen at the junctions. At the low frequencies 
which turn out to be relevant for the system dynam- 
ics, cj <C (Re[ZT]Cj)"\ we can take Re[ZT] = Re[Zo]. 
We further assume that the embedding circuit provides 
a low, real (Ohmic) impedance, Zq = R <^ Rq, such as 
would be generated by a transmission line J^ Our descrip- 
tion can easily be extended to take into account a finite 
impedance in series with the gate volt age, ^^ but since we 
take the limit Cg <C Cj this has a much weaker influence 
and so we neglect it here. 

The voltage dependence of the charging energy leads to 
resonances where the eigenvalues of the charging Hamil- 
tonian, Hch [Eq. [2], become degenerate. The charging 
energies of the states |0,/c) and |l,A:+p + 1/2), where 
p = 0,1,2,..., become degenerate at particular val- 
ues of the drain-source voltage, V = Vres , given by 
(2j9 + l)eV}i'J = Ec{l - 2ng). The energy levels near 
the p = and p = 1 resonances are illustrated in Fig. [2j 

Close to degeneracies in the charging energy the 
Josephson coupling becomes important and the pairs of 
states |0, /c) and |l,A:+p + 1/2) become strongly mixed 
forming doublets. The interaction with the environment 
can then cause the system to decay into the neighbor- 
ing doublets with lower energy. Taken together the co- 
herent evolution and decay form a cascade in which k 
increases systematically and hence a dc current flows. 
Away from the resonances the coherent evolution is sup- 
pressed and decay processes cannot take the system to 
ever higher k values so the current is also suppressed. 
Thus, the degeneracies in the charging energy lead to 
resonances in the cwcienti^^^^^ The one exception 
to this picture arises for the p = resonance where one 
can see from Fig. [2fb) that for voltages above resonance, 
V > Ec\l — 2ng\/e^ the system can move indefinitely to 



larger values of k via incoherent decay processes alone 
and hence in this case the resonance becomes strongly 
broadened on one side^i^S 



III. MASTER EQUATIONS 

Having seen how and where the Cooper-pair reso- 
nances arise, we now proceed to obtain a detailed quan- 
titative description of the charge dynamics of the SSET 
that includes the dissipation and decoherence induced 
by the electromagnetic environment. As a first step we 
use a unitary transformation method to derive an ef- 
fective Hamiltonian which provides a systematic way of 
accounting for the coherent effect of the Josephson cou- 
pling between resonant states. We then proceed to derive 
the master equation for the SSET tracing out the envi- 
ronment using the Born-Markov approximations. The 
resulting master equation can then be used to derive a 
much simpler equation that describes just the SSET is- 
land charge. 



A. Effective Hamiltonian 

Close to the p-th Cooper-pair resonance, and pro- 
vided that the voltage, V, is not too small, the eigen- 
values of the SSET charging Hamiltonian are grouped 
into doublets, {|0, /c) , |1, /c + p + V^)}' with the spacing 
between members of a given doublet much less than the 
spacing between the doublets. The main effect of the 
Josephson Hamiltonian will be to introduce couplings be- 
tween states within each doublet. Since the full system 
Hamiltonian cannot be diagonalized exactly, we treat the 
Josephson coupling as a perturbation and use a unitary 
transformation to derive an effective Hamiltonian which 
takes into account the mixing it induces between states 
within a doublet. ^^ 

We seek a unitary transformation, /7, such that, H^ = 
UHsU^ is block diagonal in the space of the doublets. 
This transformation is found as a perturbation series in 
J, and we keep only the leading order contributions (de- 
tails are given in the appendix). This results in an effec- 
tive system Hamiltonian which is block-diagonal in the 
pairs of nearly-resonant states {|0,/c),|l,/c+p + V^)}- 
Each block takes the form, 



Hi = 



E-AE 



E + AE 



(6) 



where E = —2eVk is the average charging energy of the 
resonant states, and Jp is the high order coupling between 
the states. For p > 1, 



Jp = i-if 



J" 



(2eVi^J 



2p 



(p'-Y 



(7) 



where q = 2p + 1. For the case p = 0, Jp = J. The 
splitting between resonant states, AE, is given by, 



AE = 



Ec{l - 2ng) - qeV J^ 



2q 



eVi'J q' 



1 



(8) 



The first term is the electrostatic energy difference be- 
tween the states and the second term is a correc- 
tion which arises at second order in the perturbation 
calculation^^. For p = the second order correction is 
given by JV(eF/e?). 

The block Hamiltonians are diagonalized by a rotation, 
Ua = e~*^^", where ay is the usual Pauli matrix, to give 
the eigenstates of the doublets, 

\a,k) =cosa\0,k) + sina |1, A: + (2p+i)/2) , (9a) 

\b,k) = -sina|0,/c)+cosa|l,/c + (2p+i)/2), (9b) 



and the corresponding eigenenergies, Ea^k 
Eij^k = E -^ AE\ where a is defined by 



sin 2a - 
cos 2a 



-J. 



p 



AE^' 
AE 



AE'' 
and the energy level splitting is 



AE' = sgn{AE)JAE^ + J2, 



E - AE', 

(10a) 
(10b) 

(11) 



which changes sign at the resonance. 

Note that the description of the system in terms of 
doublets is only valid within a region around each reso- 
nance. The spacing between energy levels in the doublet 
should be much smaller than the spacing between the 
doublets, \AE'\ <C eV. 



B. Born-Markov description 

We now use the block-diagonal form of the Hamilto- 
nian to derive the master equation for the SSET. We 
assume that the interaction between the SSET and the 
bath is weak,^^ R <C Rq, and that the bath has a suffi- 
ciently dense spectrum of levels that the standard Born 
and Markov approximations can be made.^^ 

Written in terms of the eigenstates of the system 
Hamiltonian, the Born-Markov master equation for the 
components of the SSET density operator, a, takes the 
formJi 



dau 



dt 






(12) 



where the tilde denotes the interaction picture and the 
sum is over only the secular^^ terms for which cj^^^ = 
^/iV, with Ujj^iy the frequency difference between eigen- 
states /i and u. The coupling tensor ^^v^i'v' is given by 



dr 



l{r) 5uv' 



/_^^fin^nfi'^ '"'" 



^fifi'^u'u^ ^'^ 




Eh' h' J^nu'^ 



. h' h' f^'^^uu'^ 



, (13) 



where k' = UkW and ^(r) = 4{e/h)^{SV{T)SV{0)) is the 
correlation function of the electromagnetic environment 
whose properties were specified in Eq. (J5]). 

The transformed operator k' (given explicitly in the 
appendix) takes the form of a power series in J, 



k' = fc(o) + A:(i) + fc(2) + 



(14) 



order term k[^/^,klJ, 



and we proceed by expanding the terms of the form 
in ^ up to second order in J. The zeroth 

is diagonal in the charge state ba- 
sis and generates dephasing of the charge states. For 
states within the same doublet, this leads to dissipative 
transitions between the eigenstates (intra-doublet transi- 
tions). The next non-zero contribution comes from terms 

of the form k[^\ klJ, , which link states from a given dou- 
blet with states in the nearest neighbor and next-nearest 



neighbor doublets, leading to inter-doublet transitions. 
The same generic description applies to all the resonances 
with p > 1, but for the p = case the states in a given 
doublet only couple to one other doublet leading to a 
slightly different form for the master equation as we dis- 
cuss below. 



To calculate the inter-doublet terms we note that close 
to resonance the inter-doublet transitions occur on a 
much larger energy scale than the spacing between levels 
in the doublet, peV ^ AE' . This allows us to sim- 
plify the calculation by ignoring the Jp terms in the 
block-diagonal Hamiltonian [Eq. (J6])] and treat the charge 
states |n, k) as the eigenstates of the system. Thus, us- 
ing the charge state basis, the inter-doublet contributions 



take the form 



dt 



da 



1/C+V2,1/C' + V2 



dt 



da^ 



0fe,l/c' + V2 



dt 



— rA/c<^Ofc,0/c' , 

(15a) 

= -{^L + Ti^ + rA/e)^lfe+V2,l/c' + V2' 
inter 

(15b) 

Ta/c 1 CrOfe,l/c' + V25 

(15c) 



Fl+T^ 



where we have assumed kuT <€. peV and the dephasing 
rate is given by 



Afc 



Sy jO) 



Afc2, 



(16) 



where uOah = —^ha = —2/\E'/h. 

Finahy, we transform the inter-doublet contributions 
[Eq. (p!5|) ] into the eigenstate basis and combine them 
with the intra-doublet terms [Eq. (p!9|) ] to obtain the fuh 
master equations of the system. For sufficiently weak 
dissipation, ^ l{r) <^ \^ah\i it is possible to simplify the 
master equation significantly by applying a further ro- 
tating wave approximation (RWA). After this approxi- 
mation the master equation takes the form, 






\^ab + ^aa 



pP+l 
^ ab 



pP+1 
^ aa 



-Jab -^^Ak)o-ak,ak', 



(21a) 



with Ak = ki — k2 for ^ifei,j7c2- The transition rates at 
the center of the resonance^^ are given by 



J 



Syjujp) 



^peVr^iJ 



h? 



J 



^4(p+l)eK 



(p) 



Sy {ujp j^l) 



(17) 



(18) 



dcr})k. 



bk' 



dt 



n 



~ -"- ab^aki,ak[ 
i/c2 ,ak'2 



- bb^bki,bk[ 



+ J- a6 ^ak2,ak' + J- 55 (^bk2,bk' + 7a6 ^ak,ak' 



K, 



- 66 



•76a 



" ^Ak)o'bk,bk' 1 

(21b) 



with cjp = 2peVris /^' For the p = case there is 
only one decay channel linking different doublets, and 
the associated energy difference is 2eV. Thus, provided 
2eV :^ /cbT, AE^ the same set of inter-doublet terms is 
obtained, but with Tl = 0. 

To calculate the terms in the master equation describ- 
ing the intra-doublet transitions we need to use the full 
eigenstates of the system [Eq. (|9|)] and take account of 
the effects of the Jp terms in the Hamiltonian. The en- 
ergy differences between the states within a doublet can 
be much smaller than those between the doublets, and so 
in this case we include the effects of a finite temperature. 
This leads to the intra-doublet equations. 



^^ak,ak 



dt 

dcfhk.bk 



dt 

dcFak.bk 



dt 



'- JbaO'bk,bk — labO'ak,ak, 

"Jab^akyttk ^ba^bk^bki 

'lab +76a 

^ak.bk ") 



(19a) 
(19b) 
(19c) 



where 7^^ gives the transition rate between the states 
|i, /c) and | j, k) within the k-ih doublet. These have the 
form 



ddr 



k,bk' 



dt 



rp _L. rP 

^ aa ^ -L 66 



^aki,bk[ 



pp+1 I pP+1 
^ aa ^ -*- 66 



0'ak2,bk' 



^L^^R , ^ A/c 



7coh + Ta/c d-ak,bk', (21c) 



where ki = k—p^ /c2 = k—p—1 and T^j are the transition 
rates between the states \i^k-\-p) -^ \j^k)^ which are 
given by 



aa 



= r^j^ = Tl cos a sin 


a, 


(22a) 


= rP+^ ^r^cos^asin^a, 


(22b) 


Tl = TLCos^a, 




(22c) 


rC'=Tncos^a, 




(22d) 


rl, = rLsm^a, 




(22e) 


Tit' =rnsm* a. 




(22f) 



lij 



2p- 



2 • 2 Sv{(^ij) 
cos asm a — 



h^ 



(20) 



We have also defined 









/coh ~ \ {2p±l 



■ 2 2 Sv(0) 

sm a cos a ^ ' 






Ak = 0, 



I 



(23) 
(24) 



For the RWA to be valid we require that incoher- 
ent decay rates are much smaher than the Bohr fre- 
quency associated with the doublets, which, for p > 1, 
results in the condition Tl <C 2\Jp\/h (since Tl is the 
largest decay rate). From Eqs. (O and (p!7|) , we find 
that, hTL/\Jp\ (X {R/RQ){eVri^J/J)^P-\ By tuning the 
gate voltage, VreJ can take on any value in the range 
< {2p-\- l)eK-es < Ecj thus, for a given value of R, the 
requirement that the RWA is valid puts a limit (which 
becomes stricter as p increases) on the maximum volt- 
age that can be considered. Our interest here is in the 
regime where the SSET charge dynamics is largely coher- 
ent, leading to sharp resonances in the current, and hence 
we naturally focus on the regime where the RWA is valid. 
Since these conditions can only be met in practicei^i^iiS^ 
for the lower values of p, we will concentrate on the 
p = 0, 1 resonances. 

The master equations bear a strong resemblance to 
those which describe the radiative cascade of quantum 
optics. ^^ In the radiative cascade a laser field drives a 
two-level atom and when the field is treated quantum 
mechanically the eigenstates of the system are atom- 
field hybrids (dressed states). Decay processes lead to 
a cascade in which photons are emitted and the laser 
state moves progressively towards lower photon numbers. 
Analogously in the SSET the states {|a, /c), |6, /c)} are 
like the atom-laser dressed states, with the island charge 
states playing the role of the atom and k like the state 
of the laser. In this case decay processes generate a cas- 
cade in which the system evolves towards states of ever 
increasing k. 



where p^ = {pla^ Pbb^ Pab)^ - Within the RWA the form 
of the matrix M, is given by 



M = 



-(r, + r^) n 



^ a 





-(r. + r?) 








where F^ = TAk=q and the other rates are given by, 
_ (rL + rfl)(l + 2cos^asin"a) , , , ^ 

^ coh ~ o "*" '*'coh "I" ^ ' 



(27) 

(28) 
(29) 

q- 
(30) 



By taking q = we obtain the equation which describes 
the evolution of the island charge. The resulting mas- 
ter equation describes a two-level system (TLS) which is 
both driven and damped. This simple two-level descrip- 
tion is all that is required to calculate the main charge 
noise properties of the SSET. 

The full density operator does not have a steady state, 
the system cascades to increasing values of k as charge 
tunnels though the transistor. The reduced equations 
derived above, however, do have a well defined steady 
state. All of the reduced coherences where q ^ are 
zero in the steady state, {p^f )ss = 0. The q = case 



aa{bb)/ 



^b{a)l(Xa + r5), and 



has the steady state: {p^ 

{Pab)ss ~ ^- Note that here and in what follows we drop 

the superscripts on p and F for the case q = 0. 



C. Effective two- level system 

Although the full master equation for the SSET is 
rather complicated, it is possible to derive a much simpler 
set by tracing over the charges that have passed through 
the SSET. Defining a set of reduced coherences, ^^ 



/5?iW = E<^'^ + 9kWli.fc). 



(25) 



and carrying out the trace over /c, we obtain a much 
simpler matrix equation. 



p^(t) = (2igeF + M)p^(t), 



(26) 



IV. AVERAGE CURRENT 

As a first application of the master equation (|2T]) , we 
calculate the steady state average current, (/)ss5 through 
the transistor. The current is determined by the rate of 
change of the number of Cooper-pairs which have crossed 
the device. 



(/),, =2e^=2eTr[fc'a]. 



(31) 



We calculate only the dominant term which comes from 
the lowest order part of the k' expansion, /c*^^^ = /c, as the 
next lowest order contribution (from /c^^^) vanishes and 
we neglect higher order contributions. The only non-zero 
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FIG. 3: (Color online) Current-voltage characteristics of the 
SSET in the vicinity of (a) the p = resonance (b) the p = 1 
resonance. The current is scaled by /o = e{rL + I'r) in each 
case. The values of the other parameters are given in the 
text. In (a) we also show in the dashed, red line the current 
calculated using the full voltage dependence of the decay rate 
Tr and the Hamiltonian. 



terms after performing the trace come from the dissipa- 
tive parts of the diagonal master equations [Eqs. (j21ap - 
(|21bp ] which give rise to the incoherent transitions. Since 
these terms only depend on the reduced coherences of the 
system, the current reaches a well-defined steady-state in 
the long time limit, 

(/)ss = 2e {{paa)ss siu^ Q^ + {Pbb)ss COS^ Ot) 

x(prL + (p + i)r^). (32) 

Within the regime where the RWA is valid, this matches 
the results obtained previousl y^^i^^i^^ using a rate equa- 
tion approach and Fermi's golden rule. 

The current near the p = and p = 1 resonances is 
shown in Fig. [3] for the typical parameter values -'^^ ^^^^^^ 
Ec = 4.Ej = 100 /ieV, T = 30 mK, a resistance for the 
embedding circuit R = 50Q and we have set rig = 0.1 
[we use these parameters throughout for the numerical 
calculations]. The full lines in Fig. [3] show the current 
calculated with the decay rates, Jp and the second order 
correction to AE given by their on-resonance values. ^^^^^ 
As the p = resonance is very broad (in comparison to 
the p = 1 resonance), we have also calculated the current 
including the full voltage dependence of the relevant de- 
cay rate and the Hamiltonian [dashed curve in Fig.[3l^a)]. 

The two resonances show rather different characteris- 



tics. The current around the p = resonance is broad 
and highly asymmetric. This is because in this case 
purely dissipative processes can generate a dc current 
for V > Vres (incoherent Cooper-pair tunneling^^) as can 
be seen from the energy level diagram in Fig. [2jb). The 
current for the p = 1 resonance is much closer to the stan- 
dard Lorentzian form of a resonance, the small amount 
of asymmetry still present arises from the intra-doublet 
transitions. For V < Kes, relaxation between the lev- 
els of the doublets (controlled by Jabiba)) hinders current 
flow whilst for V > Vres it helps it. This leads to a small 
asymmetry in the current as a function of voltage which 
is only removed when the temperature is sufficiently high 
such that Jab — Iba and a Lorentzian shape is recovered. 
Extending our calculation to include the regime where 
the RWA is no longer valid, we find that the current peaks 
at the resonances become suppressed. This is because 
outside of the RWA the system is unable to build up the 
coherence between charge states necessary for current to 
flow. This is consistent with Ref. 28 where this effect was 
studied in detail. 



V. CHARGE NOISE SPECTRUM 

We now turn to consider the noise properties of the 
SSET close to the Cooper-pair resonances. The charge 
noise spectrum provides a fingerprint of the subtle inter- 
play of coherent dynamics and dissipation in the system 
whilst also controlling the back-action that the SSET will 
exert when it is used as a measuring device. 

The charge noise spectrum is given by. 






{dn'{t)5n'{0))e"^'dt, 



(33) 



where Sn' 



' — {n') is the fluctuating part of the trans- 
formed island charge operator and the averages are taken 
over the steady state. The operator n' can be written as 
a series expansion in terms of powers of J (as we did for 
k')^ n' = n^^^ + n^^^ + . . ., as described in the appendix. 
The dominant contribution to the spectrum comes from 
the zeroth order terms in the expansion, n^^^ = n, and 
so we consider these terms first. The higher order terms 
in n' give rise to weaker features in the spectrum which 
we go on to calculate in Sec. IV CI 



A. Triplet structure 

The charge noise spectrum is readily obtained using 
the quantum regression theorem^^ which allows us to cal- 
culate the behavior of (n(t)n(O)) for t > by using the 
equation of motion for {n{t)) with appropriately modified 
initial conditions. The starting point for the calculation 
is the effective master equation for the SSET charge state 
[Eq. ([261)] with q = 0. 




■2 crq 



3 ^^ 
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i^ 



hjj/ Jp 




huj/ Jp 




FIG. 4: (Color online) Charge noise spectrum as a function 
of frequency and voltage around the (a) p = and (b) p — 
1 resonances. The resonances occur at the point where the 
central peak disappears. 



We begin by defining a series of projection operators 

P^J=I]b■,^)(^,^^ (34) 

k 

with the property {pij) = (pij) ss which then allow us to 
calculate the fluctuating part Spij = pij — (p^)ss- The 
correlation function for the charge written in terms of the 
projectors is, 

{6n{t)Sn{0)) = ^ nijnim{Spij{t)Spim{0)) , (35) 

i,j,l,'m 

where the matrix elements are given by the transfor- 
mation of n into the eigenstate basis, Uaa = sin^ a, 
7255 = cos^ a and Uab = nija = sin a cos a. 

Using the regression theorem and the initial conditions 
for the projectors, 

{Spij{0)5pim{0)) = 6im{pji)ss - {Pij)ss{plm)ss^ (36) 

we obtain the spectrum, 

p 

^ ~^ i pop 



2^2 / \P55/ssrcoh 



(Paa/ss-L coh 



{co-co,,)2 + rl^ (^_^,j2 + r2^^ 



(37) 



where Fpop = Ta -\- T^. Examples of the spectrum as 
a function of frequency and voltage around the p = 



FIG. 5: (Color online) Charge noise, Snn, as a function of 
A^ for (a) Q = 100 MHz and (b) Q = 4 GHz oscillator. SoUd 
(blue) lines are Snn{^), dashed (red) lines are Snn{—^)- The 
inset in (a) shows the region around the minimum where the 
difference between the two curves is most pronounced. 



and p = 1 resonances are shown in Fig. SI In both 
cases the spectrum consists of three Lorentzians. This 
triplet structure is exactly what is expected for a coher- 
ent TLS in the presence of dissipation^*^ The central 
peak around zero frequency arises due to incoherent tran- 
sitions between eigenstates, its height is determined by 
\^aa — ^66 1 7 the difference in average charge between the 
I a) and \b) eigenstates. This peak disappears when the 
system is tuned to resonance since the eigenstates are 
equal mixtures of charge states. 

The sidepeaks arise at cj = ±uJab due to coherent os- 
cillations between the eigenstates. The heights of the 
sidepeaks are controlled by the steady state populations 
of the eigenstates, which is what provides the asymme- 
try between the negative and positive frequency peaks, 
and the width is given by the rate at which the coher- 
ences decay, Fcoh [see Eq. ([30]) ]. Note that the spacing 
of the sidepeaks is much larger for p = 0, the splitting 
Jp=o ^ 25 GHz, compared to Jp=i ^ 1 GHz for p = 1. 
This means that finite temperature effects are not as im- 
portant in the p = resonance. 



B. Back-action 

We now consider the effect that the charge noise spec- 
trum has on another system coupled to the SSET island. 
For concreteness we consider the case of a resonator with 



a charge-position coupling, 

He = Xn{c^c^), 



(38) 



where c is the resonator lowering operator, and A the 
coupling strength. For weak coupling, and provided that 
the resonator does not start to self-oscillate, ^^^^^ linear 
response theory can be used to calculate the effect of the 
SSET on the resonator dynamics J^i^S This back-action of 
the SSET on a resonator can be useful, providing in some 
cases a way to cool the resonato r^^i^^ whilst also provid- 
ing a way of measuring the asymmetry in the SSET's 
noise properties. ^^ 

Within the linear response regime, the SSET acts on 
the resonator like an additional thermal bath and its ef- 
fect can be characterized by an effective damping, jba^ 
and an effective thermal occupation number, uba- For a 
weakly damped resonator of frequency 1], these are given 
byii^iO 



jBAm = X^[Snnm-Snn{-m. 



riBAi^) 






(39) 
(40) 



Measurement of jba and uba allows both frequency 
components of the charge noise spectrum, Snn{^^)^ to 
be inferred. We will explore the back-action effect of the 
SSET on resonators in two different regimes of frequency, 
Q = 100 MHz and ft = 4 GHz, corresponding to different 
physical realizations of the resonator. The 100 MHz fre- 
quency is typical of a nanomechanical resonator^^ whilst 
the 4 GHz frequency matches that of superconducting 
striplines^ or LC resonators.— We focus on the p = 1 
resonance from now on since the sidepeaks of the p = 
resonance are too wide to have any significant effect on 
even a high frequency LC oscillator. 

Figure [5] shows the charge noise as a function of the 
detuning from resonance for Q = ±100 MHz and Q = 
±4 GHz. At 100 MHz the spectrum is very symmetric, 
Snn{^) ^ Snn{—^) bccausc the sidepeaks have very little 
weight at this frequency. In this case maximum asymme- 
try is achieved at the center of the resonance {AE = 0) 
when the spacing of the outer peaks is minimized. In 
contrast, at 4 GHz the spectrum is highly asymmetric as 
here the sidepeaks cross through this frequency. 

The curves for Snn{^^) are not simply reflections of 
each other, as would be expected for a classically driven 
TLS^° or other resonances in the SSET J^ii^ This asym- 
metry occurs for the same reason that the current peak is 
not a simple Lorentzian; the intra-doublet decay rates are 
not symmetric between ±A^. Similar effects are seen in 
a driven TLS when^^ the temperature dependence of the 
relaxation rate is taken into account. The behavior of 
the driven TLS where Snn{^,^E) = Snn{-^,-^E) is 
recovered in the high temperature limit where the intra- 
doublet transition is saturated, jab = Iba- We can un- 
derstand how the intra-doublet decays lead to this asym- 
metric behavior in terms of effective temperatures. The 



intra-doublet decays drive the eigenstate populations to- 
wards an equilibrium distribution which corresponds to 
the temperature of the bath, this is constant and al- 
ways positive. However, the inter-doublet rates drive the 
SSET to an equilibrium point whose effective tempera- 
ture varies strongly with AE and, in particular, changes 
sign when AE = 0. The competition between these two 
behaviors causes the asymmetry in this system. ^^^^^ 

In Fig. [6] we plot the damping, jba^ for Q = 100 MHz 
and Q = 4 GHz. The small asymmetry in the low fre- 
quency noise spectrum gives the damping a very small 
magnitude, but for the high frequency resonator it is 
much larger. The lack of symmetry in the noise spec- 
trum, Snni^^^E) ^ Snn{—^^ —^E)^ Icads to quitc dif- 
ferent magnitudes for the damping at ±AE, with the 
anti-damping peak suppressed by the intra-doublet de- 
cays for both the high and low frequency cases. 

When a resonator is coupled to the SSET its steady 
state is determined by a combination of the back-action 
of the SSET and the influence of the rest of the res- 
onator's surroundings which are in thermal equilibrium 
at a temperature T, and give rise to a damping rate 7ext- 
The overall occupation number of the resonator, ntot, is 
given by the average , ^^^^^ — 



^tot 



7ext^ + iBAriBA 
Text + 7BA 



(41) 



where n = {ex.-p{hQ/kBT) — 1)~^. Thus at a given T, the 
SSET can be used to cool the resonator provided pro- 
vided Uba < n. Such cooling is important for nanome- 
chanical resonators with frequencies in the MHz range 
as even at temperatures T below 100 mK they will still 
contain a large number of thermal quanta^i^i^i^ 

Figures ^c) and [6l^d) show the behavior of uba (plot- 
ted as I/uba to emphasize the behavior near mini- 
mas) around the resonance at 100 MHz and 4 GHz. For 
ft = 100 MHz, the minimum in uba occurs when the sys- 
tem is tuned directly to resonance where the central peak 
in the noise spectrum vanishes. For the typical device 
parameters we have chosen (given in Sec. HV|) . we find 
a minimum of uba ^ 37 which corresponds through the 
expression n^A = {e^'p{hQ/kBTBA) — ^)~^ to an effective 
temperature Tba ^ 29 mK, only slightly lower than the 
bath temperature, T = 30 mK. This result is, however, 
dependent on T through the intra-doublet rates and the 
relative cooling potential does improve at higher bath 
temperatures, for example Tba ^ 37 mK for T = 50 mK 
and Tba ^ 50 mK for T = 80 mK. 

The main problem with using the SSET tuned to a 
Cooper-pair resonance to cool a mechanical resonator is 
the spacing of the peaks in the noise spectrum. For ef- 
fective cooling, we need the frequency of the resonator 
to match the separation of the peaks. This cannot be 
achieved at low frequencies in this device since the min- 
imum intra-doublet spacing (and hence minimum peak 
splitting) 2Jp/his in the GHz range. Trying to engineer 
a device where this splitting was much smaller would 
lead to a deterioration in the effectiveness of the cooling 
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FIG. 6: Back-action damping, ^ba [(a)- (b)] and effective occupation number, uba [(c)- (d)] of an oscillator weakly coupled to 
the SSET. (a) and (c) are for O = 100 MHz, (b) and (d) for O = 4GHz. 



because of the effect of thermal noise on the SSET it- 
self: the asymmetry in the peaks would be reduced and 
they would also be broadened (this would take the sys- 
tem outside both the regime where the RWA is valid and 
also the resolved sideband limit where optimal cooling 
can be achieved^^^^^). In contrast, for a driven TLS^^ 
the potential for cooling is much greater as it is possible 
to tune the drive (which corresponds to the Josephson 
coupHng in our system) and it is the difference between 
this frequency and the level separation of the TLS(both 
of which can be ^ k^T /K) that sets the spacing of the 
sidepeaks in the corresponding noise spectrum. 



C. Higher order spectral features 

The full charge noise spectrum also contains contribu- 
tions at other frequencies far from cj = 0. These arise 
from the higher order contributions in the perturbative 
expansion of n' . In this section we investigate the most 
significant of these high order terms, which arises from 
the {n^^\t)n^^\{f)) term in the expansion. The form of 
n^^^ is calculated in the appendix. 



.(1) 



E 



J 



2peV 



\0,k){l,k+y2\ 



J 



2{p + l)eV 



\0,k){l,k-y2\+h.c., (42) 



for p > 1, in the case p = the first term is not 
present. When this is transformed to the eigenstate ba- 



sis it contains terms proportional to \i, k) {j, k + q\, with 
q = ±p, ±{p+l), forp = only terms q = ±1 are present. 
As an example we calculate the corresponding spectrum, 



S. 



.(i)„(i) 



/CO 
(n(^)(t)n(i)(0))e'"*dt, 
-CX) 



(43) 



for the p = 1 resonance (note that in the steady state 
(n^^^) = and so we do not need to work with the opera- 
tor Sn^^''). To do this we use the regression theorem and 
the reduced master equations [Eq. (|26|) ] with the values 
of(7 = ±l,±2. 

We find that the spectrum consists of triplets of peaks 
centered on the frequencies hu = 2qeV with the side- 
peaks separated by the intra-doublet level spacing in each 
case. The triplets with g = ±1 and q = ±2 simply differ 
by a constant prefactor and a slight modification to the 
decay rates and so we concentrate here on just the q = ±1 
case. The spectrum around huj = ±2eV is shown in Figs. 
[Tfa) and[7fb). The peaks in this part of the spectrum 
arise from first order Josephson coupling between states 
at the relevant frequency differences. We show the rele- 
vant transitions and their frequencies in Fig. [Tl^c). Since 
this spectrum occurs at higher order in the perturbation 
theory the magnitude of 6'^(i)^(i) is much smaller than 
it was for the zeroth order spectrum (Fig. H]). The cen- 
tral peak in both Figs.[71^a) and[71^b) is the same because 
the transitions which give rise to these peaks always link 
corresponding states within the two doublets (i.e. an a 
state with an a state or a 6 state with a b) for both the 
positive and negative frequency processes as shown in 
Fig. [Tfc). However, the sidepeaks of the two triplets are 
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FIG. 7: (Color online) 5'^(i)^(i) (cj) in the vicinity of the p = 1 resonance (a) and (b) show the spectrum at frequencies around 
hw — zL2eV respectively. A schematic illustration of the processes giving rise to the spectral features is shown in (c). This 
illustrates the inter-doublet transitions and the corresponding frequencies of the features they give rise to. 



quite different. The weights of the sidepeaks are propor- 
tional to |n^^.p(Pii)ss where n^^- = ^j^ {i, k ^ q\ n^^^ |j, k) 
and i{j) is the initial (final) state for the relevant tran- 
sition. This means that each of the four sidepeaks has a 
different combination of matrix element and population. 

Far from the resonance only one of the three peaks is 
present in each triplet. In this region the eigenstates are 
very close to the charge states and since the Josephson 
effect only couples the states |0, k) and |1, /c + 1/2) we find 
that each state is only coupled to one other. Very close 
to the resonance all of the peaks appear, the eigenstates 
are mixtures of charge states and so transitions between 
all of the states in the doublets can occur. 

The features seen in this part of the spectrum arise 
because the system is not as simple as a true TLS, they 
arise from couplings between different doublets and hence 
require more than two energy levels. The frequencies 
at which the features in part of the spectrum occur, 
~ 100 GHz, are much larger than the range that can be 
probed with a stripHne resonator. However, it might be 
possible to observe the noise at this frequency in a dif- 
ferent kind of experiment in which the SSET is instead 
coupled to another mesoscopic conductor such as a SIS 
junctionals 



VI. CONCLUSIONS 

We have analyzed the quantum dynamics of the SSET 
tuned close to Cooper-pair resonances. An effective 
Hamiltonian for the SSET was derived, this exploits 
the separation of the energy levels into doublets and 



accounts for the Josephson coupling between resonant 
states. We then derived the master equations for the 
system including the effect of the electromagnetic envi- 
ronment using the Born-Markov approximations. We cal- 
culated the current in the vicinity of the resonances and 
find, in accord with previous studies, strong peaks at the 
resonances. 



Calculating the charge noise spectrum in the vicinity of 
a resonance, we found that the spectrum is dominated by 
a triplet of peaks centered on zero frequency: a structure 
typical of a driven, damped, TLS. However, the detailed 
form of the triplet in this case differed in important re- 
spects from a standard classically driven TLS because of 
the intra-doublet transitions. An experimentally realiz- 
able method of measuring the quantum noise spectrum 
is to couple the SSET to a resonator and measure the 
back-action. Using a linear response approach we found 
that the effects of the intra-doublet relaxation can be 
observed in the asymmetry of the damping rate, 7b A, 
and effective occupation, uba^ of the resonator. It would 
also be possible to cool a low frequency oscillator using 
a Cooper-pair resonance, though not to the ground state 
as the large separation of the sidepeaks in the noise spec- 
trum, 2Jp/h, strongly limits the minimum occupation 
number that could be achieved. We also carried out a 
detailed calculation of the higher order triplet features 
that arise in the charge noise spectrum. These features 
occur because of the Josephson coupling between differ- 
ent doublets. 
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Appendix A: Transformations 

1. Transforming the Hamiltonian 

We begin with the system Hamiltonian Hs = i^ch + 
Hj as defined in Eqs. (|2]) and (|3]). We find a unitary 
transformation, H' = UHW (we drop the subscript S 
in this appendix), such that H' only contains diagonal 
elements and those which couple resonant states. To do 
this we define the projector onto the k-th doublet as 

Pk = |0, k) (0, fe| + |1, fc + p + 1/2) (1, ^ + P + V2| , (Al) 

then the condition on H' above becomes PkH'Pk" = for 
k ^ /c'^ We then write the transformation a.^ U = e**^ 
with S = S^ and treat Hj as a perturbation. This allows 
us to write a series expansion for S in terms of J, 



S = 5(°) + J5(i) + J25(2) + . . . 



(A2) 



There are infinitely many transformations which satisfy 
the condition on H' and so to uniquely specify U we 
choose that S should not have matrix elements within 
the doublet, PkSPk = 0. We note that ^^^^ = 0, since 
for the case J = we need no transformation. This then 
allows us to find the transformed form of the Hamiltonian 



as 






which we can then write as a power series 
where we have grouped terms by order in J; 



1 r 
2! 



iJS^^\[iJS^^\H^^] 



(A3) 

(A4) 

(A5a) 
(A5b) 

(A5c) 



These can then be used to construct the expression for S 
term by term and so build up the effective Hamiltonian. 



As an example we give the calculation of S'^^^. We begin 
by noting that PkH'^^^Pk" = for all n and k 7^ k" which 
gives, 

= PkHjPk'' + Pk{iJS^'^H,i, - H,i,iJS^'^)Pk», (A6) 

from which we can calculate the matrix elements of iS'*^^-' 
as 



{0,fc|iS'(^)|l,A;- 1/2} 



J 



2peV 
J 



2ip+l)eV 



= Gi, (A7) 

-G2, (A8) 



and the obvious conjugates. Gi and G2 then become 
the natural small parameters of the perturbation theory. 
These expressions can then be used to calculate S'^'^^ and 
so on. This allows the transformed Hamiltonian to be 
found up to any order, for example we find that the sec- 
ond order corrections to the diagonal elements are given 

by 



(0,A;|ff(2)|o,A;) 



J2 2q 
'eVq^ 



1 



(A9) 



We also need to be able to calculate the off-diagonal cou- 
pling between resonant states, Jp. This appears to or- 
der 2p -\- 1 in the perturbation expansion. The method 
outlined above becomes very cumbersome at high orders 
and so to calculate Jp for higher orders we introduce the 
level-shift operator —i2^ 



R{z) = Hj^Hj- 



-Hj 



Hch 

+ Hj-^-Hj-^Hj + ..., (AlO) 

where Qk = 1 — Pk- Which then allows us to calculate 
Jp = (0, A: I R{E) |1, A; + p + 1/2) • This gives the expres- 
sion found in Eq. ([7|). However, we have checked explic- 
itly that this approach gives the same expression for Jp^i 
(the highest order used in the main text) as that obtained 
using the perturbation expansion. 



2. The transformed operators 

We also need to transform the operators k and n. To 
illustrate the method we give an explicit calculation of 
k' . This is done in exact analogy with the calculation of 
H\ We first write k^ as a power series. 



fc' = fc(0) + fc(l) + fc(2) + . . _ . 



(All) 



where. 





fC — fV, 


(A12a) 




fc(i) = [iJ5«,fc], 


(A12b) 


fc(2) = 


[ij25(2)^fc] + i.[^J5(l)j^J5(l)^fc]" ^ 


(A12c) 
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which then allow us to calculate k'-^^ as 
fcW = ^^Gi|0,fc){l,fc+V2| 

k 

+ G2|0,A:)(1,A:- V2I + /1.C, (A13) 



and /c^^^ is given by 



^(2) ^ {2p + l)G,G, ^ ^1^^ ^^ ^^^ ^ ^ ^1 _ |g^ ^^ ^^^ ^ _ ^1 ^ ^^^^^ 

+ (G? + G2)^(|l,fc + V2)(l,fc + V2|-|0,fc)(0,fc|). (A14) 

I 



A similar calculation can be performed to calculate the 
power series for n' = n^^^ + n'^^^ + n^^^ + . . ., and we find 
n^ = n, 

n^i) = ^ Gi |0, k) (1, /c + V2I-G2 |0, fe) (1, /c - V2|+/^.c., 

^151 



n^^^ = G1G2 Y^ (|0, fe) (0, A: + 1| - |0, A:) (0, fe - 1| + /i.c) 

+ (G?+Gl)^(|l,A: + V2)(l,fc + V2|-|0,fc)(0,A:|). 

(A16) 
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